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Abstract. The three-point correlation function of cosmic shear, the weak distortion of the images of distant galaxies 
by the gravitational field of the inhomogeneous matter distribution in the Universe, is studied here. Previous work 
on three-point statistics of cosmic shear has mainly concentrated on the convergence, or on aperture measures of 
the shear. However, as has become clear recently for the two-point statistics of cosmic shear, the basic quantity 
that should be used is the correlation function: first, it is much easier to measure from observational data, since 
it is immune against complicated geometries of data fields (which contain gaps and holes, e.g. due to masking); 
second, all other (linear) two-point statistics can be expressed as integrals over the correlation function. The 
situation is the same for the three-point statistics. However, in contrast to the two-point correlation function, 
the invariants (with respect to rotations) of the shear three-point correlation function have not been employed 
yet. Here we consider the transformation properties of the shear three-point correlation function under rotations. 
We show that there are four complex linear combinations of components of the three-point correlation function, 
which we shall call 'natural components', since they are multiplied just by a phase factor for arbitrary rotations, 
but do not mix. In particular, their moduli are invariant under rotations and thus (non-linear) invariants of the 
three-point correlation function. In terms of these natural components, the invariance of the statistical properties 
of the shear field under parity transformations are easily obtained. Our results do not apply only to cosmic shear, 
but also to other quantities with the same mathematical properties - that of a polar. For example, practically 
every relation derived here applies also to the polarization of the cosmic microwave background radiation. 

Key words, cosmology: theory - gravitational lensing - large-scale structure of the Universe - cosmic microwave 
background 

1. Introduction 

The weak gravitational lensing effect by the large-scale matter distribution of the Universe, called cosmic shear, has 
long been recognized as a unique tool to study the statistical properties of the cosmological (dark) matter distribution, 
without referring to luminous tracers of this distribution (Blandford et al. 1991; Miralda-Escude 1991, Kaiser 1992, 
1998, Jain & Seljak 1997, Bernardeau et al. 1997, Schneider et al. 1998, van Waerbeke et al. 1999; Bartclmann & 
Schneider 1999; Jain et al. 2000, White & Hu 2000; see Mellier 1999 and Bartelmann & Schneider 2001 for recent 
reviews). Owing to the smallness of the effect, its actual measurement has only fairly recently been achieved, nearly 
simultaneously by several groups (Bacon et al. 2000; Kaiser et al. 2000; van Waerbeke et al. 2000; Wittman et al. 
2000). This breakthrough became possible due to the usage of wide-field optical cameras and the development of 
special-purpose image analysis software specifically designed to measure the shape of very faint galaxies and to correct 
their ellipticity for effects of PSF smearing and anisotropy. By now, several additional cosmic shear measurements 
have been reported (Maoli et al. 2001; van Waerbeke et al. 2001; Rhodes et al. 2001; Bacon et al. 2002; Refregier 
et al. 2002; Hammerle et al. 2002; Hoekstra et al. 2002), both from the ground and from HST imaging, partly with 
appreciably larger sky area than the original discovery papers. 

In all of these papers, the cosmic shear signal detected was one related to the two-point correlation function of the 
shear, or some function of it, such as the shear dispersion or the aperture mass. To measure higher-order statistical 
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properties of the cosmic shear, the quantity of data must be larger than for the second-order measures. It has been 
pointed out by a number of authors (e.g., Bernardeau et al. 1997; Jain & Seljak 1997; Schneider et al. 1998; van 
Waerbeke et al. 1999; Hamana et al. 2002) that the third-order statistics (e.g. the skewness) contains very valuable 
cosmological information, such as the density parameter Q m . In particular, the near-degeneracy between erg and fl m 
in two-point cosmic shear statistics (see, e.g., van Waerbeke et al. 2002) can be broken if the three-point statistics is 
employed. Encouragingly, Bernardeau et al. (2002a) have reported the detection of a third-order statistical signal in 
their cosmic shear survey. 

Apart from the larger difficulty to obtain a measurement of the third-order statistics of the cosmic shear, there is 
also the problem of an appropriate statistical estimator for the third-order shear. Whereas for the second-order, the 
statistically independent shear measures are known, we are not in this position for the third-order shear statistics. We 
shall briefly summarize the situation for the two-point statistics, and explain why the three-point shear statistics is 
substantially more complicated in Sect. 2 below. In Sect. 3 we shall then define the components of the shear three-point 
correlation function (3PCF) and study their transformation behavior under spatial rotations. From that, we shall then 
find in Sect. 4 the natural components of the shear 3PCF, which can be considered analogous to the natural components 
of the two-point correlation function of the shear. These components are 'natural' in the sense that they have the 
simplest behavior under rotation transformations: each of the four complex natural components is just multiplied 
with a phase factor when an arbitrary rotation is applied, which in particular means that the moduli of these natural 
components are invariants under rotations. We shall discuss the importance of these natural components in Sect. 5, 
where we also outline the perspectives of future work that can be based on the use of these natural components. In an 
appendix, we shall consider the projection of the shear onto several particular reference points, defined by the various 
centers of a triangle. 

We want to point out that all the relations derived in this papers are not confined only to cosmic shear. In fact, 
this paper investigates the three-point correlation function of a polar - a polar is a two-component quantity which 
transforms under a rotation of the coordinate frame by a phase factor e 2lip . Alternatively, a polar can be viewed as the 
trace-free part of a symmetric 2x2 matrix. In the case of cosmic shear, this matrix is the Hessian of the deflection 
potential. Another polar of great cosmological importance is the polarization of the cosmic microwave background. 
Hence, all the results presented below do equally well apply to the three-point correlation function of the CMB 
polarization. 

2. Motivation 

For the two-point cosmic shear statistics, the basic quantities are the two-point correlation functions. Given a pair 
of points, Xi, and the Cartesian components of the shear 7 M (Xi) there (i = 1,2, fj, = 1,2), one projects the shear 
along the direction ip connecting these two points by defining the tangential and cross component, 7t and j x by 
7ti + hxi = -[7i (Xj) + i7 2 (Xj)] c~ 2lv . Then, one forms the correlation functions £«(#) = (7ti7t2) and £ X x(#) = 
(7x17x2}, where the average is an ensemble average over all pairs of points with separation 0. Even more useful are 
the linear combinations £±(0) — £tt(6*) ±£ x ><(#). Another combination which one may be tempted to take is (7ti7x2), 
but this changes sign under parity transformations and thus should vanish. All other two-point statistical measures of 
the cosmic shear, such as the shear dispersion in a circle or the aperture mass dispersion, can be expressed as integrals 
over these two correlation functions (e.g., Crittenden et al. 2002; Schneider et al. 2002a). From a practical point of 
view, the determination of the shear correlation functions is also most convenient, as they can be measured also in 
data fields of complicated geometry (as normally data fields are, due to masking). The two correlation functions £± 
can be expressed readily in terms of the power spectrum of the mass distribution in the Universe. 

Compared with this situation, a proper measure of the three-point shear statistics is much more difficult to define. 
First we note that the two-point function £ + = (77*} (where '*' denotes complex conjugation) can be defined without 
any reference direction; this is not the case for any three-point function of the shear, since with three two-component 
quantities alone, no tri-lincar scalar can be formed. Hence, one needs to project the shear components. In contrast to 
the case of two points, where there is a unique choice of the reference direction, this is no longer true for a triangle: 
there is not a single 'natural direction' defined in a triangle; in fact, there are several of those (see the Appendix). 
Therefore, it is not a priory clear how to define 'useful' components of the three-point correlation function. One might 
ask, for example, whether there are similar 'invariant' combinations of the components of the shear 3PCF as there are 
for the two-point function (£±). 

Given these difficulties, it is not surprising that the work on the three-point statistics of cosmic shear has been 
relatively sparse. Bernardeau et al. (1997) and van Waerbeke et al. (1999) consider the 3PCF of the surface mass 
density as reconstructed from the shear measurements. Whereas possible in principle, the fact that real data sets 
have gaps and holes makes the reconstructed mass map susceptible to systcmatics due to the geometry. Schneider et 
al. (1998) suggested to use the aperture mass (Schneider 1996) as a cosmic shear statistics for which the third-order 
moment is readily calculated directly from the shear data. However, as is the case for the shear dispersion, one needs to 
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Fig. 1. Definitions of the geometry of a triangle. The ipi Fig. 2. The orthocenter of a triangle, determined by the 

are the orientations of the sides of the triangle relative to intersection of the three altitudes. The angle between the 

the positive xi-direction, the <pi are the interior angles of altitude of Xi and the side X3 is easily obtained from the 

the triangle, which are related to the ipi by ([!]), provided right triangle formed by this altitude and the two sides Xi 

the orientation of the three points is as displayed here. and X3. 



cover the data field with (circular) apertures which presents again a problem in case of gaps in the data. Recognizing 
this, Bernardeau et al. (2002b) defined a particular component (or, more precisely, a particular linear combination 
of components) of the 3PCF that is readily measured from observational data, calculated its expectation value from 
numerical ray-tracing simulations of Jain et al. (2000) and successfully applied it to the VIRMOS-DESCART survey 
in Bernardeau et al. (2002a). 

In addition to the component of the three-point shear correlator considered by Bernardeau (2002b), the other 
components (there are a total of 8) may contain equally valuable information about the bispectrum of the mass 
distribution. In addition, it is easily seen that all measures of the third-order shear statistics can be expressed as 
integrals over the shear 3PCF. They are easiest to determine from real data and shall therefore be considered as 
the basic quantities. In this paper we will derive the 'natural' components of the shear 3PCF, by considering their 
transformation behavior under rotations. 



3. The shear three-point correlation function, and the centers of a triangle 

3.1. Definition of the 3PCF 

Consider three points X;, 1 < I < 3, and define their difference vectors xi = X3 — X2, X2 = Xi — X3, X3 = X2 — Xi, 
so that xi + X2 + X3 = (see Fig. [[]). Each of the three difference vectors will be written as x; = (xi cos ipi, xi sin ipi), 
so that if 1 is the orientation of the i-th side of the triangle. Furthermore, we define <pi as being the interior angle of 
the triangle at the corner X;. In order to have general relations, in this paper we will use oriented angles, i.e. we 
will attach to each angle a sign indicating its orientation. More precisely, we will define <j>\ G (— 7r, +ir) to have the 
same sign of the cross product X2 x X3 (where axb:= 0162 — 0261), and similarly for 02 and 03. Note that, since 
xi x X2 = X2 x X3 = X3 x xi (which follows from the vanishing of the sum of the x; , or in a more geometric way, from 
the fact that each of these cross products equals twice the area of the triangle), all angles 0/ will be either positive or 
negative; in particular, they all will be positive if the closed path from Xi to X2 to X3 to Xi goes around the triangle 
counter-clockwise. We also observe that this convention for the angles implies that the sum of the internal angles of 
the triangle will be ±7r depending on the orientation; however, this ambiguity will not generally play a role, since all 
relations for angles are defined modulo 27r. In the following we will call a triangle positively (respectively, negatively) 
oriented if the sum of its internal angle is +ir (—it). The relation between the 0/ and the tpi is given by 

¥>3 - ¥>2 = 7T - 01 ; fl ~ <^3 = 7T - 02 ; <p2 ~ <Pl = 7T - 03 ■ (1) 

Let 7 M (X/) be the Cartesian components of the shear at point X;; they are defined in terms of the deflection 
potential V»(X) by the differential operation 7l = [d 2 ip/(dX 2 ) - d 2 tp/ (dX$)]/2, 72 = d 2 ij)/{dXi dX 2 ). We define the 
Cartesian components of the shear 3PCF as 



7awa(X1,X2,X 3 ) = (7 m (Xi)7^(X2)7a(X3)} , 



(2) 
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where we have made use of the fact that the shear field is assumed to be a homogeneous random field, so that the 
3PCF is invariant under translations; therefore, 7 m „a does depend only on the separation vectors xj. Although we also 
assume that the random field is isotropic, the Cartesian components 7 m; ,a of the 3PCF do not depend just on the 
X[ = |x;| (see below). 

For any reference direction Q, we can define the tangential and cross components of the shear, jn and 7; x , 
respectively, at point X; relative to this direction by 

7 (X / ,0) = 7t(X ; ,0) +i 7 x(X i ,6) := -7(X)e- 2iC ' - - [71 (X,) + i 72 (X,)] e" 2 ^ , (3) 
or explicitly in terms of the components, 

7„(X I ,C0 = -JW2C07,(X I ), (4) 

where R^ v {ip) = cos V5 + e M „ sin tp is the rotation matrix, 6^ v the Kronecker delta, and £12 = 1 = —£21, en = = £22; 
here and in the following, we identify the tangential component with fj, = 1, and the cross component with fi = 2 
when we write transformation equations like (ji|), and we employ the Einstein summation convention. We notice that, 
because of the rotation behavior of the shear, we only need to specify the projection direction Q modulo n. 

If the reference directions Q are defined in terms of the position vectors of the vertices of the triangle, or the side 
vectors, and thus rotate in the same way as the triangle as whole, then the tangential and cross components of the 
shear are invariant under rotations and translations of the triangle. For example, the direction Q could be chosen as 
the direction ipi of the opposite side of the triangle. The 3PCF of these projected shear components will then depend 
only on the |x/|, 



^,1(X1,X2,X 3 ) = -i? Ma (2Cl)i? J , / 3(2C2)-RA 7 (2C3)7a/3 7 (xi,X2,X3) . (5) 



The orientation of the triangle becomes important if the arguments of the 3PCF are written as the three side lengths; 
in this case, the triangle is defined only up to a parity transformation. Hence, in all relations involving only the side 
lengths of a triangle, we will assume that the triangle has positive parity. We shall see in Sect. 4 how the shear 3PCF 
behaves under parity transformations. 

In Eq. (j5|), the directions Q are arbitrary. For a given set of three points, i.e., for a given triangle, there are several 
natural choices for the reference directions; we shall discuss those in the next subsection. One choice was already 
mentioned above, namely the direction ipi of the side opposite to the corner X;. We shall label the corresponding 
3PCF with the superscript 's' (for 'side'), 

%v\( X l> X *i X *) = -^"( 2( ^l) Rufj{^2) i?A 7 (2^3)7a/3 7 ( x l> x 2,X3) ■ (6) 



3.2. The centers of a triangle, and the corresponding shear projections 

In addition to the side projection mentioned above, there are several other natural choices for the directions along 
which the shear can be projected. For each triangle, one can define a number of 'centers'; the direction of the vector 
connecting the point X/ with one of these centers can be used to define convenient components of the 3PCF. The four 
most important centers are: (1) The centroid of a triangle; it is the point where the side-bisectors intersect; (2) the 
incenter (center of the incircle), which is the intersection of the three angle-bisectors of the interior angles 4>i; (3) the 
circumcenter (center of the circumcircle) , which is the point of intersection of the three midperpendiculars, and (4) 
the orthocenter, which is the intersection point of the altitudes. 

We now define the tangential and cross components of the shear relative to the direction of the line connecting 
the point X/ with one of these centers. Since the line connecting X/ with the orthocenter (the point 'H' in Fig.||) is 
perpendicular to the side vector x;, we find for the projection of the shear relative to the direction of the orthocenter 
(labeled with superscript 'o') simply by setting Q = tpi + ir/2 in (|4|), 

ll o) = ~li s) ■ (7) 
Therefore, we obtain the shear 3PCF for this projection from Eq. (|5|) 

T^Ja^i, £2,2:3) = -^l x (x 1 ,x 2 ,x 3 ) . (8) 

Next we consider the center of the incircle. As this is given by the intersection of the angle-bisectors, we obtain for 
the direction Ci of the line connecting Xi with the incenter (see Fig.0) £i = ip 3 + cj>\/2; hence, 



7, 



( in )(X 1 ) = -i?^(2^ 3 + ( /» 1 )7 ! ,(X 1 ) 
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Fig. 3. The incenter of a triangle is given by the intersec- Fig. 4. The centroid of a triangle is the intersection point 
tion of the three interior angle-bisectors. of the three side-bisectors. 



Using Eqs. (|) and (|), one finds that the relation between the Cartesian components of the shear and those defined 
with respect to the ortho-center is 



7 M (X / )=i?^(-2 W )7^° ) (X ; ) , 
so that 

7 ( in) (Xi) = -R^ptpa - 2<p 1 + ^ih^Xi) • 

Since 2tp3 — 2ipi + <t>\ = —2ir + 202 +<f>i = - 
7r, one finally obtains 



(in) 



R b 



h - h)Rvf3 (03 - 0l)i?A 7 (01 - 4>2)T a 



(o) 
07 



(9) 



b 3 , where we used Eq. (|l]) and the fact that the sum of the 0/ is 



(10) 



where the result for the other two points were obtained from the one at Xi by cyclic permutation of the indices. It 
should be noted that in this equation, as in all the following, the angles 0; are functions of the sidelengths Xi (and 
the orientation), and not independent geometrical quantities; hence, an equation like ( |l0| ) contains only the x\ as 
independent variables. 

Next we consider the projection onto the center of the circumcircle. If d denotes the direction of the line connecting 
Xi with this center, then the shear components with respect to this center (denoted by the superscript 'out') read 



^^(2Ci-2^ 1 ) 7 (° ) (X 1 ) 



(11) 



7 (°"t) = _i^(2Ci) 7 ,(Xi) 

where in the second step we used Eq. (||). From Thales' theorem (see Fig. |]) one finds that Ci — ¥?3 = ±tt/2 — 03, where 
the sign ± depends on the orientation of the triangle. Hence, £i — <p\ = ±7r/2 — 03 + (p^ — (pi — = F 7r /2 + 4>2 ~ 03, and 
we see that the sign ambiguity does not play a role because of the rotation properties of the shear [cf. comment after 
Eq. (0)]. For the other two points, the corresponding relations are obtained by cyclic permutations, so that 



7%f = R„a{2h - 203)^(203 - 200^(20! 



20 2 )7, (O) 



a/37 



(12) 



Finally, we consider the projection on the center of mass of the triangle, or the centroid (denoted by the superscript 
'cen'), which is the intersection of the three side-bisectors. Denoting by ^3 the angle between the side bisector of side 
X3 and X3 , we obtain from Fig. |] that 

. , xi . , X1X2 . , 

sin V> 3 = 7— sin 02 = 7 sin 3 , 

h 3 h 3 x 3 

where we made use of the sine-theorem, and 

1 



2^2 



is the length of the side-bisector of side x 3 . From the same figure we also obtain that 
cos tp 3 



hl + xl/4 



x\ 



x\ 
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h 3 x 3 



2h 3 x 3 
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Fig. 5. The circumcenter of a triangle is given by the in- 
tersection of the three midperpendiculars. According to 
Thales' theorem, the side X2 subtends the angle 202 as 
seen from the circumcenter; from that, the angle between 
the lines connecting the points X; with the circumcenter 
are easily obtained by noting that all three of them have 
equal length. 



Fig. 6. The centers of the escribed circles are intersections 
between one interior angle bisector and two exterior angle 
bisectors. These three escribed circles are tangent to one 
side of the triangle and the extensions of the two other 
sides. 



so that 



{x\ - x\ ) 2 - A:x\x\ sin 2 3 . {x\ - x\)xiX2 sin 4> 3 



cos ^ = 5 sm 2 ^ = ■ (13) 

Analogous relations are obtained for the other two points by cyclic permutations of the indices. The shear components 
projected toward the centroid are then 

7 ( con )(X0 = i?^(2^)7 (s) = -iV(2^)7 (o) , (14) 

where the components of the rotation matrix can be directly obtained from Eqs. (|l3|). For completeness, we give the 
3PCF for the centroid, 

7^ = -iW 2 ^i)^M2V>2)iW 2 ^3hS 7 ■ (15) 



4. Invariant components of the shear three-point correlation function 

The normal way to choose 'good' components of a multi-index quantity like 7 m „a is to look for its behavior under 
coordinate transformations. Recall the situation for the shear two-point correlation function: the correlator (7^7^) (x) 
of the Cartesian components of the shear contains two terms, one which is independent on the phase ip of x, the other 
behaving as cos(4<y9). The coefficients of these two terms are £+ and respectively, which are the invariants, and 
therefore the natural choice for the components of the two-point correlation function of the shear. 

Similarly, we can consider the behavior of the 3PCF under rotations. Given that no linear scalar can be built from 
the Cartesian components of the 3PCF alone we cannot expect to find a combination of components which is invariant 
under general rotations. However, as we shall see, there are linear combinations of the components of the 3PCF which 
have a simple behavior under rotations, namely multiplication by a phase factor. Owing to this simple transformation 
behavior, we shall term them 'natural' components. 
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Let now 7 m „a be the components of the 3PCF of the shear measured with respect to one particular choice of 
direction(s); choosing different projection directions, which differ from the old ones by Q, the 3PCF becomes 

= (7;(XiK(x 2 )7i(x 3 )) = iW2Ci) iM2&) Rx-yiXsha^ ; (16) 

note that this expression differs from Eq. <S) by a minus sign since in Eq. (|^) the transition from Cartesian to tangential 
and cross components was included as well. To investigate more closely the action of the rotation operator in Eq. (pf), 
we transform this equation into one which appears more familiar: We dchnc the eight-component quantity 

T := (7ttt,7ttx,7txt,7txx,7xtt,7xtx,7xxt,7xxx) 



and analogously V for the transformed components; then, the rotation described in Eq. (16) can be written as T' = RT, 
where R is a 8 x 8 matrix. The components of R are triple products of trigonometric functions. Since it describes a 
rotation, R is unitary and one expects that the eigenvalues of R have an absolute value of unity. They can in fact be 
obtained as 

Ag = exp (±2i[Ci + C 2 + Cs]) ; A$ = exp (±2i[-Ci + C2 + Cs]) ; 

\fl = exp (±2i[d - C2 + Cs]) ; Ag = exp (±2i[d + C2 - Cs]) ■ (17) 

Note that this result is not very surprising: either the eigenvalues are ±1, or they have to appear in the above form, 
i.e., they have to occur as pairs of complex conjugate numbers with absolute value of unity. The dependence of the 
eigenvalues on the rotation angles is in fact a natural one. 



4.1. The natural components of the 3PCF 

We see that under the general rotation described by Eq. ([l6|), the eigenvalue +1 does not occur; in other words, there 
is no (linear) combination of the components of the 3PCF that is invariant under the transformation (|l6|). However, 
for some special rotations, the eigenvalue +1 occurs. As a first example, we consider the case Ci +C2 + C3 = 0; as shown 
in Sect. 3.2, this case actually is encountered for transformations of the shear components between the orthocenter, 
the incenter and the circumcenter. Then, the two eigenvalues A^ are +1, and the corresponding eigenvectors can be 
found to be 

E{ 0) = (1,0, 0,-1, 0,-1, -1,0); E$ ] = (0,1, 1,0, 1,0, 0,-1). (18) 

Hence, under rotations of this kind, we expect that the two combinations 

r^ 0) = 7ut - 7txx - 7xtx - 7xxt ; r^ 0) = 7 ttx + 7 txt + 7 xtt - 7 XXX (19) 

of the components of the 3PCF are both invariant. It will turn out to be very useful to combine these two invariants 
into a single complex quantity, — + ir 2 ' ■ This quantity, however, can also be written in a different form, 
namely 

r«« = (7(X 1 ) 7 (X 2 )7(X 3 )) , (20) 

where we again consider 7 = 7t + 17 x as a complex quantity. Written in this form, it is obvious that a rotation with 
Ci + C2 + Cs = keeps r(°) invariant. 

Next we consider the case that — Ci + C2 + C3 = 0- Such a rotation occurs if one transforms the shear components 
from the incenter (or the out- or orthocenter) to the center of one of the escribed circles (see Fig.^ for an explanation); 
in this case, A ^ = 1, and from the corresponding eigenvectors one can again construct two combinations Ij , T^ 1 "*, 
of the components of the 3PCF that stay invariant under such a rotation, and they can as well be combined into a 
single complex quantity. The same procedure can then be repeated for rotations with Ci — (2 + C3 = an( l those with 
Ci + C2 — Cs — 0; these correspond, e.g., to the transformation from one of the aforementioned centers to the centers 
of the other two escribed circles. The corresponding invariants read: 

= 7ttt ~7txx +7xtx +7xxt +i[7ttx +7txt - 7xtt +7xxx] = (7*(Xi)7(X 2 )7(X 3 )) ; 
T (2) = 7ttt + 7txx - 7xtx + 7xxt + i [7ttx - 7txt + 7xtt + 7xxx] = (7(^1)7* (X 2 )7(X 3 )} ; (21) 
r (3) = Tttt + 7txx +7xtx - 7xxt + i [-Tttx +7txt +7xtt + 7xxx] = (7(Xi)7(X 2 )7*(X 3 )} . 

Now, each of these r( Q ) is invariant only under special rotations, as the derivation above has shown. But the remarkable 
finding here is that, under a general rotation, the different r' Q ) do not mix, but they are just multiplied by a phase 
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factor; indeed, since the rotation ( [Tq ) implies that 7(Xi) — > 7'(Xi) = 7(Xi)e 2l ^ 1 , and similarly for the other two 
points, we see from the final expression in Eqs. ( |20| ) and (|2l]) that the transformed invariants read 

(r<°>)' = exp(-2i[Ci +C2 + C 3 ])r^ = A 2 0) r(°) , and (I»)' = \ { 2 a) , (22) 

where the A2 have been defined in Eq. (|l7|). The transformation (|22| ) is indeed remarkable, though in hindsight not 
all that surprising. The fact that the T^'only transform amongst themselves justifies that they are called natural 
components of the three-point correlation function of the shear. Note in particular that Eq. $2!^) implies that the four 
absolute values |r( Q )| are independent under general rotations; they are therefore the (non-linear) invariants of the 
shear 3PCF. 

Of course, if desired, the original components of the shear 3PCF can be obtained from the natural components; 
the inverse of Eq. (^lj) reads: 

7ttt = (rW+r[ J ) + r[ 3 » + r<°')/4; 7 txx = f-r« + rf> + if> -if>) /4 ; 



7xtx = ( r« - if > + - r<°>) /4 ; 7x xt = (if 5 + r?> - if - t®) /a 

7tt> 



(r« + if) - if > + r<°>) /4 ; 7txt = (r« - r< 2 > + r 2 3) + r 2 °>) /4 ; (23) 

7x« = (-if ) + if> + r 2 3) + if) /4 ; 7x x x = (r« + r 2 2) + if > - r 2 °>) /4 . 

4.2. Parity transformation 

A parity transformation is obtained by mirroring a given set of three points with respect to any straight line. Such a 
transformation has the following effects: if in the original triangle the three points Xi, X2, X3 are ordered counter- 
clockwise with respect to, say, the center of the incircle, in the transformed triangle they will be ordered clockwise. 
As a result, under a parity transformation we need to replace every angle with its opposite. 

In order to understand the behavior of the 3PCF under parity, let us suppose that the triangle is flipped along 
a line perpendicular to the side X3 (since the 3PCF is invariant upon rotation, the direction used for the flipping is 
irrelevant). We first observe that this mirror symmetry is equivalent to interchanging the points Xi and X2, or the 
sides x\ and xi- Second, such a flipping of orientation keeps the tangential component of the shear invariant (if the 
direction relative to which the shear components are measured are as well subject to the mirror transformation - like 
it happens when shear components are defined relative to one of the centers of the triangle), but changes the sign 
of the cross-component; hence, this transformation implies P7 = 7*. Together, these two effects therefore imply that 
under a parity transformation, 

P [7/u* (zi.afe, a*)] = n 7 ^ A (x 2 ,xi, x 3 ) , where II = (-1)^+ A+1 , (24) 

the parity, is +1 if all of the indices of 7 are t's, or two x's occur, otherwise it is negative. Components of the 3PCF 
with II = +1 are called even, those of negative parity odd components. The action of the parity transformation is 



most easily expressed in terms of the natural components of the 3PCF; as can be seen from Eq. (21), the real part of 
each of the r'"-* is composed of even components, the imaginary part of odd components. Therefore, 

P [r<°5(a;i, o 2 , a*)] = (r (0) ) * (n, o 3 , x 2 ) = (r<°>) * (a*, x lt x 3 ) = (r (0) ) * (a*, x 2 , Xl ) , 

p [r^zr, aM)] = (r^y ( Xl ,x 3 ,x 2 ) = (r^y (x 2 ,xx,x 3 ) = (r^y (x 3 ,x 2 , Xl ) , 

p [r( 2 ) ( Xl , x 2 , x 3 )] = (r^y ( X2 ,x u x 3 ) = (r^y (x 3 ,x 2 , x x ) = (r< 3 ))' (o 1; o 3 ,o 2 ) , (25) 

P [r< 3 >(a:i, x 2 , x 3 )] = (r^* (x 3 ,x 2 ,x 1 ) = (r^)* {x u x 3 ,x 2 ) = (r^* (x 2 ,x 1 ,x 3 ) . 

Each of these parity transformations contains a permutation of the arguments with negative signature. One can also 
consider cyclic permutations of the arguments; since 

l fiV \(Xl,X2,X 3 ) = r fuXn{x2,X 3 , X l) = i yx l Mu(x 3 ,Xl,X 2 ) , (26) 

one obtains for the natural components 

rW{ Xl , X2 ,x 3 ) = rW( X 2, X3 , X i) = rW{ X3 , Xl , X2 ) , 

r (1) (oi,o 2 ,o 3 ) = r (3) (02,03,01) = r (2) (o 3 ,oi,o 2 ) , 

r^(0i,02,0 3 ) = r«(0 2 ,0 3 ,0i) = r^(0 3 ,0i,02) , (27) 

r (3) (oi,o 2 ,o 3 ) = r (2) (o 2 , o 3 ,oi) = r (1) (o 3 , oi,o 2 ) . 
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Hence, r^ -* is invariant under cyclic permutations of the arguments, whereas the other three natural components of 
the shear 3PCF transform into each other under such permutations. The relations (|2^) imply that only one of the 
three functions T 1 --^ {x\, ie 2 , £3), 1 < fc < 3, is independent, the other two can be obtained by cyclic permutations of 
the arguments. 

The parity transformations have an immediate consequence for triangles where two sides are equal, say x\ — X2', 
namely, from Eq. (|24|) one finds that 

7ttx (^1,^1,^3) = = 7xxx(zi,xi,x 3 ) ; ltxt[xi,x u x^) = ~ j xtt (x 1 ,x 1 ,x 3 ) . (28) 

This implies that for x\ = X2, T^ ^ and F( 3 ^ have no imaginary part, and those of and I^ 2 ) have equal magnitude 
but opposite sign. Furthermore, for equilateral triangles, all odd components of the shear 3PCF vanish, in which case 
the natural components become purely real. 



4.3. Two-point correlation function revisited 

The foregoing formalism can of course also be applied to the two-point correlation function. In that case, the eigenvalues 
are A^ ^ = e ±2l ^ 1=t ^ 2 ), and the invariant combinations are = 7t t — 7xx + i (7tx + 7xt) and = 7a + 7xx + 
i(7tx — 7xt). The two-point correlator is, however, special in the sense that after a parity transformation, the two 
points can be brought back into the old positions with a rotation. This then implies that the imaginary components 
of the r^") vanish if they are measured in the only reference frame that makes sense for two points - namely, the line 
connecting them; and so T^ ) = r' 1 ) = £ + . Note that if different projection directions are taken, the imaginary 
parts of the r( Q ) are not zero. 



4.4. Generalization 



The discussion on obtaining the natural components of the shear 3PCF immediately suggests how to generalize it to 
the natural components of higher-order correlation functions. Here, we shall give the results for the four-point function 
(also see the recent work by Takada & Jain 2002a for a detailed consideration of the kurtosis of the cosmic shear field) ; 
further generalizations are straightforward to obtain: 

If "ffiuXn is the four-point correlation function of the shear (for notational simplicity, we skip the arguments of this 
function), and we consider a general rotation of the directions relative to which tangential and cross components are 
defined, the transformation of the shear four-point correlation function reads 

t£„a„ = (^(X 1 )7i(X a ) 7 UX3K(X 4 )) = iW2Ci)iM2Ca)i2A-r(2Ca) R^Wlc^s ■ (29) 
In analogy to the treatment for the 3PCF, there are now eight complex conjugate pairs of eigenvalues, which are 



exp(±2i[Ci + C2 + C 3 + C4]) 



A 



,(0) 
n,2 

= exp(±2i[Ci-C2 + C3-C 4 ]) 

fl = CX P (±2i[Cl+C2+C3-C4]) 
,(6) 



1,2 

(3) 
1.2 

(5) 
1,2 

(7) 
1,2 



exp(±2i[Ci + C2-C3-C4]) ; 
exp(±2i[-Ci + C2 + C3-C 4 ]) 
exp(±2i[Ci + C2-C3 + C4]) ; 
exp(±2i[-Ci + C2 + C3 + C4]) 



(30) 



v M - exp(±2i[Ci-C2 + C3 + C4]) 
The corresponding natural components are then 
r(°) = ( 7 (X 1 )7(X 2 )7(X 3 ) 7 (X 4 )} ; 
r« = <7(X 1 ) 7 *(X 2 ) 7 (X 3 ) 7 *(X 4 )} ; 

= < 7 (X 1 ) 7 (X 2 ) 7 (X 3 ) 7 *(X 4 )) ; 
r( 6 ) = ( 7 (X 1 ) 7 *(X 2 )7(X 3 ) 7 (X 4 )} ; 
For reference, we shall write down the first of these explicitly, 

r'°' = 7 tttt + 7 xxxx — 7ttxx — 7txtx — 7txxt — 7xttx — 7xtxt — 7xxtt 

+ i [Ttttx + 7ttxt + 7txtt + 7xttt — 7tx x x — 7xtx x — 7x xtx — 7x x xt] 
These natural components of the shear four-point correlation function will transform under the rotation 



r (i) 

r (3) 
r (5) 
r (7) 



( 7 (X 1 ) 7 (X 2 )7*(X 3 ) 7 *(X 4 )} 
( 7 *(X 1 ) 7 (X 2 ) 7 (X 3 ) 7 *(X 4 )) 
( 7 (X 1 ) 7 (X 2 ) 7 *(X 3 ) 7 (X 4 )} ; 
( 7 *(X 1 ) 7 (X 2 ) 7 (X 3 ) 7 (X 4 )) . 



(31) 



(32) 



like 



(rw)' 



A (a) p(a) 



(33) 



These relations make it obvious how generalizations to higher-order correlations can be obtained. 
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4.5. The estimator of Bernardeau et al. (2002) 

In their paper, Bernardeau et al. (2002b) considered, for the first time, a specific shear 3PCF, which they then applied 
successfully to observational cosmic shear data (Bernardeau et al. 2002a). Here, we shall write the Bernardeau et al. 
(2002b) estimator in our notation. They considered one side of a triangle (say, X3) as the reference direction, and 

project the shear at all three points along this direction. Hence, their projected components (here written as 7tj> <) 
read 

7 ( 3) (X ; ) = -iV(2^ 3 )7,(X;) = R^3 - 2 W )7^ S) (X ; ) , (34) 

and, of course, 7 (3) (X 3 ) = 7 (s) (X 3 ). The 3PCF defined by Bernardeau et al. (2002b) then becomes 

Pjf) = ([ 7t (3) (X 1 )7t (3) (X 2 ) + 7 W(X 1 ) 7 (f ) (X 2 )] 7^ 3) (X 3 )) . (35) 

Since 

7t (3) (X 1 ) 7t (3) (X 2 )+ 7 ( 3) (X 1 ) 7 ( 3) (X 2 ) 

= [R lv (2tpa - 2<p 1 )R lx {2v 3 - 2ip 2 ) + R 2v {2^ - 2tp 1 )R ax (2<p 3 - 2y 2 )] 7 ( s) (Xi) 7 ^ s) (X 2 ) 
= i?, A (2^i-2^ 2 ) 7 ( s )(X 1 ) 7 i s) (X 2 ) , 

one obtains 

V™ = R^fa)^ , (36) 
provided the triangle formed by the points X/ has the orientation defined at the beginning of Sect. 3. 

5. Discussion and outlook 

The natural components of the shear 3PCF provide a generalization of the corresponding natural components £± of 
the two-point correlation function. Given that upcoming cosmic shear surveys will be substantially larger than the 
current ones, it is obvious that the 3PCF will be measurable in the future with a similar accuracy as the two-point 
correlation function in current surveys; in fact, a first significant measurement of the 3PCF has been reported in 
Bernardeau et al. (2002a). Therefore, it is worth to explore the dependence of the 3PCF on various parameters, as 
has been done for the two-point function. In particular, the detailed study of the interrelations between £± and the 
underlying power spectrum, as presented in Crittenden et al. (2002) and Schneider et al. (2002a), shall be generalized 
to the 3PCF. Of course, this will be substantially more difficult from a technical point of view, and will be deferred 
to future work. Nevertheless, we can outline a few aspects of what can be expected from such work. 

In close analogy to the two-point correlation function, one can expect that the natural components can be more 
easily calculated from the bispectrum of the underlying mass distribution than the individual components of the 3PCF. 
As is the case for £±, one can expect that T^ ' probes the bispectrum in a different way than the T^ k \ 1 < k < 3. In 
particular, the different natural components will have a different dependence on cosmological parameters. 

The natural components of the 3PCF are not independent of each other; provided that the shear indeed is due to a 
surface mass density field, there should be integral relations which interrelate them. Again one should note the analogy 
with the two-point correlation function, where £ + can be obtained as an integral over £_ and vice versa (Crittenden 
et al. 2002; Schneider et al. 2002a). The interrelations between the components of the 3PCF provide a redundancy 
which can be profitably combined to reduce the noise in real measurements. 

If the shear is not solely due to a surface mass distribution, the shear field may contain a B-mode contribution (e.g., 
from intrinsic alignments of the galaxies from which the shear is measured) . In the case of the two-point correlation 
function, the presence of a B-mode can be probed from integral relations between the two correlation functions £± 
(Crittenden et al. 2002; Schneider et al. 2002a); it is expected that similarly in the case of the 3PCF, the integral 
relations between the natural components will be modified in the presence of a B-mode. 

All linear three-point statistics of the shear can be expressed in terms of the shear 3PCF. This is obvious from the 
fact that the bispectrum of the surface mass density can be expressed in terms of the 3PCF; on the other hand, all 
linear three-point statistics are linearly related to the bispectrum, and can therefore be expressed directly in terms of 
the 3PCF. In particular, the third-order aperture mass statistics (Schneider et al. 1998; van Waerbeke et al. 1999) can 
be expressed as an integral over the 3PCF. It remains to be seen whether the third-order aperture mass statistics is as 
useful for a separation of the shear field into E- and B-modes as it is the case for the two-point statistics (Crittenden 
et al. 2002; Schneider et al. 2002a; for applications in cosmic shear surveys, see e.g. Pen et al. 2002; Hoekstra et al. 
2002). 
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Concerning practical measurements of the 3PCF, the procedure is straightforward: from a given catalog of galaxy 
images with position vectors Oi and measured ellipticities e*, triplets are selected. For each such triplet, the sides of 
the corresponding triangle can be calculated, and for practical reasons the largest side be called £3. The three points 
are then labeled X/, 1 < I < 3 in the (unique) way such that the orientation of the points is as described at the 
beginning of Sect. 3, and the longest side connects X1 and X 2 . Since the connecting vectors X; need to be calculated 
anyway it is simplest to project the ellipticities along the sides of the triangle, or, with opposite sign, towards the 
orthocenter, without having to calculate any trigonometric function. The eight (=2x2x2) triple products of the 
projected ellipticities are calculated and summed up in eight three-dimensional bins. Those are conveniently labeled 
by a scale 2:3, and two shape parameters, q\ = xi/x^, 172 = X2/X3, so that the grid of bins runs as < (71,2 < 1, 
< X3 < 00. After summing up over all triplets of points, for each of the eight components the sums in the bins 
are divided by the number of triplets contributing, which then yields an estimate of the 3PCF. Those can then be 
combined into the natural components for further analysis. 

For a quantitative analysis of the three-point statistics of cosmic shear in relation to predictions of cosmological 
models, one can either use integrated properties of the shear 3PCF - such as the aperture mass or the integral of r( B ) 
over an elliptical region as done in Bernardeau et al. (2002b), or consider the (more noisy, but much more numerous) 
estimates of the natural components of the 3PCF directly. The first method yields one-dimensional functions of the 
three-point shear statistics depending on a scale parameter, and are thus very convenient for graphical displays. In 
contrast, using the full shear 3PCF employs multi-dimensional data which is difficult to display, but contains all 
the information about third-order statistics from the data. Hence, even though the signal-to-noise in each bin of the 
3PCF can be small, its overall information content cannot be smaller than that of any of the integrated quantities, 
and should therefore be employed in extracting cosmological information from the measurements. To obtain a reliable 
figure-of-merit for the comparison of the observational results with model predictions, one needs to know the covariance 
matrix of the shear 3PCF - this, however, is a 64-components quantity depending on 6 arguments, and will be very 
difficult to obtain analytically (see Schneider et al. 2002b for the difficulties of obtaining the covariance of the two-point 
correlation functions). Corresponding covariances of integrated quantities may be slightly easier obtainable. On the 
other hand, it is quite likely that figures-of-merit will have to rely heavily on future ray-tracing simulations through 
N-body-generated cosmological mass distributions, as in, e.g., van Waerbeke et al. (1999), Jain et al. (2000) and 
Bernardeau et al. (2002a). 

After this paper was finished, two recent preprints were posted on the Web which discussed very similar issues 
(Zaldarriaga & Scoccimarro 2002; Takada & Jain 2002b). In the latter one, the components of the correlation 
function were calculated from ray-tracing simulations; in particular they verified that all of the eight components of 
the 3PCF are non-zero in general and thus contain cosmological information. 
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